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Abstract 

We give an argument for deriving analytically the infrared "Abelian" dom- 
inance in a gauge invariant way for the Wilson loop average in SU(2) Yang- 
Mills theory. In other words, we propose a possible mechanism for realizing the 
dynamical Abelian projection in the SU(2) gauge-invariant manner without 
breaking color symmetry. This supports validity of the dual superconductivity 
picture for quark confinement. We also discuss the stability of the vacuum with 
magnetic condensation as a by-product of this result. 
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1 Introduction 



Quark confinement is still an unsolved and challenging problem in theoretical particle 
physics, as is well known. Though the approach to this problem is not unique, we have 
a promising scenario for explaining quark confinement, the so-called dual supercon- 
ductivity picture [2] for the vacuum of the non-Abelian gauge theory [1] . This scenario 
proposed long ago is intuitively quite appealing. Indeed, the relevant data supporting 
the validity of this picture have been accumulated by numerical simulations especially 
since 1990 and some of the theoretical predictions [3,4] have been confirmed by these 
investigations: infrared Abelian dominance [5], magnetic monopole dominance [6] and 
non- vanishing off-diagonal gluon mass [7] in the Maximal Abelian gauge [8], which 
are the most characteristic features for the dual superconductivity. In spite of these 
facts, the theoretical justification is not yet reached to a satisfactory level. In this 
paper, we demonstrate analytically the infrared Abelian dominance in the Wilson 
loop average for a large Wilson loop of the SU(2) Yang-Mills theory [1]. 

For this purpose, we need to answer how to define and extract "Abelian" part V M 
from the original non-Abelian gauge field which is responsible for the area decay 
law of the Wilson loop average. This must be done without spoiling gauge invariance. 
The conventional Abelian projection [3] is too naive to realize this requirement. In 
sec. 2 and 3, we achieve this goal by using a non-Abelian Stokes theorem and a 
non-linear change of variables (NLCV). 

At the same time, we must answer why the "remaining" part X M (X M = V M ) in 
the non-Abelian gauge field ^ decouple in the low-energy (or long-distance) regime. 
To answer this question, we argue in sec. 4 and 5 that the remaining part X M acquires 
the mass Mx to be decoupled in the low-energy region. The fundamental mechanism 
for the X M -mass generation is that the gauge-invariant composite operator X^ develops 

a non- vanishing vacuum expectation value (x^, in other words, the gauge-invariant 

dimension-two condensation takes place, i.e., ^X^ ^ 0. In fact, it was recently 
proposed in [9] and [10, 11] that dimension-two vacuum condensations composed of 
gluon field are relevant to the realization of quark confinement and the existence of 
mass gap in Yang-Mills theory. We also discuss some implications of dimension-two 
condensate ^X^ for the low-energy description of Yang-Mills theory. 

Moreover, in sec. 5 we point out that the existence of such a condensation stabi- 
lizes the vacuum of the Savvidy type [12] with the magnetic condensation by eliminat- 
ing a tachyon mode causing the Nielsen-Olesen instability [13]. This is a by-product 
of the above result. The stability of the magnetic vacuum is desirable for the magnetic 
monopole dominance. Thus the non-perturbative Yang-Mills vacuum is characterized 
by two vacuum condensations, i.e., the condensation (x^ ^ and the magnetic con- 
densation (H) 0, both of which realize the vacuum energy lower than that of the 
perturbative vacuum. 
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2 Non- Abelian Stokes theorem for the Wilson loop 
operator and introduction of color field 

For the non- Abelian gauge potential, srf^ix) = £/ A (x)T A , the Wilson loop operator 
Wj^^C) for a closed loop C is defined by 



W^{C) := Jf- x Xx ^expUg ^dx^^x)] 



(2-1) 



where & denotes the path-ordered product and the normalization factor jV is equal 
to the dimension of the representation R, to which the probe of the Wilson loop 
belongs, i.e., N = dim(l#) = tr(l^). Then the Wilson loop average W(C) := 
(W / £ /(C)) YM , i.e., the vacuum expectation value of the Wilson loop operator, is given 
by the functional integration: 

W(C) := (WAC)) YM =Z Y l / V^exp {iS YM )W^{C) 

_ JV^ex V (iS YM )W^C) 

JV^exp(iS YM ) ' 1 ' } 

where S Y m is the Yang-Mills action. This expression is rather formal. For this to be 
a precise definition, we must specify the gauge fixing procedure to give a well-defined 
functional integration measure T>£/ ll , as will be discussed later. 

The Wilson loop operator is rewritten into a surface-integral form which is called 
the Non- Abelian Stokes theorem (NAST). We adopt in this paper the Diakonov- 
Petrov version [14] of NAST 1 which does not include neither the path ordering along 
the loop nor the surface ordering in sharp contrast to the other versions [19] of NAST, 
at the price of an additional integration over all gauge transformations of the given 
non-Abelian background field. 

The Diakonov-Petrov NAST is nothing but a path integral representation of the 
Wilson loop operator. This representation is obtained according to the usual proce- 
dure of obtaining the path-integral representation: i) partitioning the closed loop C 
into N infinitesimal segments, ii) inserting the complete set at each partition point, 
iii) taking the limit N — > oo appropriately. As the complete set to be inserted, we use 
the coherent state which is described by introducing an auxiliary vector field n(x). 
The vector field n(x) is hereafter called color field in relation to the Yang-Mills theory 
by the reason to be clarified later. 

In what follows, we use a notation F ■ G :— F A G A = 2tr(FG), F 2 := F ■ F, and 
(F x G) A := e ABC F B G c = -2itr(T A [F, G]) with the normalization for the Hermitian 
generators T A of the Lie algebra £f of the gauge group G: tr(T A T B ) = \5 AB '. For 
the gauge group G = SU(2), the color field n(x) is the unit vector field with three 
components, i.e., 

n(x) = (ni(x),n 2 (x),n s (x)), n(x) ■ n(x) := n A (x)n A (x) = 1, (2.3) 



1 We adopt the coherent state representation for a derivation of the Diakonov-Petrov NAST as 
given in sec. Ill of [16] for SU(2) and in [17] for SU(N),N > 3. See also [18] for more information. 
It was claimed that some care must be taken in using this version of NAST [15]. 
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and the path-integral representation reads 

W*{C) = J d/2c[n] exp UgJ j d^{\x{a^U^) + ig^ti(^Ud^)}^ , (2.4) 

where J = |, 1, |, 2, ■ ■ • is the index, say "spin", characterizing the representation R, 
to which the probe of the Wilson loop belongs, and dfic[n] is the product measure of 
the normalized invariant Haar measure dfi[n(x)] on SU(2)/U(1) = S 2 at a spacetime 
point x: 

d/jLc[n] := Y\ dfi[n(x)], dfi[n(x)] = d 3 n(x)5(n(x) ■ n(x) — 1). (2.5) 

Here the unit vector field n(x) is defined through an SU(2) matrix field U(x) by 

n(x) :=n A (x)a A = W(x)a 3 U(x) (A = 1,2, 3), (2.6) 
with Pauli matrices a^(A = 1,2,3). Then, by using the Euler-angle representation 

U{x) = e ia ^ a ' 3 / 2 e if3( - x ^ CT2 ^' 2 e i ' y ^ CT ' i ^ 2 , (2.7) 
the color vector field n(x) is expressed by two Euler-angle fields a(x) and /3(x): 

n(x) = (tia(x))a=i,2,3 = (sin.{3(x) cosa(x), sin/?(x) sina(x), cos@(x)), (2.8) 
and an explicit form of the Haar measure is given by 

2 J + 1 

dfi[n(x)} = — sin P(x)dP(x)da(x). (2.9) 

4:71 

Since the argument of the exponential is Abelian, we can use the ordinary Stokes 
theorem to rewrite the line integral to the surface integral: 

WAC) = f d^ s [n] exp Ugl f dFS^oA , (2.10) 

J l Z J S:dS=C ) 

where S is an arbitrary surface spanned on the loop C, the antisymmetric tensor G^ v 
is the curvature defined by 

G^ v {x) =d lt [n(x) ■ sf v {x)\ ~ d v [n(x) ■ ^{x)\ - g~ l n{x) ■ (9 M n(z) x d v n(x)), (2.11) 

and dns[n] is the product measure over the surface S: 

d^ s [n] := Hdfi[n(x)]. (2.12) 

x&S 

See e.g., [16,37] for details of the derivation. 
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Figure 1: The relationship between the original Yang-Mills (YM) theory and the master 
Yang-Mills (M-YM) theory. The master Yang-Mills theory has a larger (local and global) 
gauge group G than the original gauge group G of the original Yang-Mills theory and 
becomes equivalent to the original Yang-Mills theory after a constraint (new MAG) is 
imposed. The resulting gauge theory is denoted by M-YM at new MAG with the gauge 
group G'. 



3 Non-linear change of variables for gluon fields 

Since the loop C can have arbitrary shape and arbitrary location in spacetime, the 
surface S spanned by the loop C may sweep the whole spacetime. In view of this, 
the color field n(x) must be introduced over all the spacetime points. This is easily 
achieved by inserting the unity, 

1 = J Vn[n] = j Vn8{n ■ n - 1) := JJ J[dn{x)]5(n{x) ■ n{x) - 1), (3.1) 

into the functional integration. Therefore, we arrive at the expression: 

W(C) =Z^ M J Vfi[n] J P^exp (iS YM )W^(C) 
fVfi[n]JV^ex V (tS YM ) W*{C) 



J Vfi[n] J exp (zSym) 
where we have introduced the reduced Wilson loop operator W^(C) defined by 

J 



(3.2) 



W*(C) = exp lig- ( d 2 S^G,A , (3.3) 

I 2 JS:dS=C J 



and the new partition function Zym has been defined by inserting the unity: 1 
/ Vfi[n] into the original Yang-Mills partition function: 



Zym = J T>n[n] J V^exp(tS Y u[^}). (3.4) 



4 



At this stage, the color field n(x) is regarded as an auxiliary field introduced into the 
Yang-Mills theory in addition to the gauge field stf^ix). We call this modified theory 
the master Yang-Mills theory, which is written in terms of ^(x) and n(x). Thus 
we regard the vacuum expectation value of the Wilson loop in Yang-Mills theory as 
that of the reduced Wilson loop in the master Yang- Mills theory. The master Yang- 
Mills theory has more independent degrees of freedom than the original Yang-Mills 
theory For a while, we put this issue aside, until we will discuss how to reduce the 
master Yang-Mills theory to the original Yang-Mills theory shortly after introducing 
the non-linear change of variables. See Fig. 1. 

We proceed to perform the (non-linear) change of variables of the original gauge 
field ^ by making use of the color field n(x). This will help us to clarify which 
variables are responsible for the area law of the Wilson loop average. Given a color 
field n(x), the Yang-Mills gauge field s^^ix) can be cast into the equivalent form: 

=0^ =(n ■ ^)n + =0^ - (n • 

= (n • ^)n + (n • n)^ — (n • £/^)n 
= (n • £$y)rt + n x x n) 

= (n • s^y)n - g^n x d^n + g~ l n x (<9 M n + gstf^ x n) 

= (n • <0^)n + g^d^n x n + g~ l n x D^[stf]n, (3.5) 

where we have used only the relation n(x) ■ n(x) = 1 in the second equality and 
introduced the covariant derivative in the last step: 

D^[£/]n(x) := d^n(x) + g^(x) x n(x). (3.6) 

Thus the Yang-Mills gauge field ^(x) is decomposed as 

^(x) = c^(x)n(x) + g^d^x) x n(x) + X M (x), (3.7) 

where we have used the identification: 

c M (x) = n(x) ■ ^(x), X M (x) = g~ l n(x) x D^]n{x). (3.8) 

The first term on the right-hand side of (3.7) is denoted by C M (x) := c M (x)n(x). C^(x) 
is parallel to n(x) and is called the restricted potential. The second term is denoted by 
B M (a;) := g~ 1 d^ i n{x) x n(x). M^(x) is perpendicular to n(x) and is called the magnetic 
potential. For later convenience, we define V M (x) by V M (x) := C M (x) + B M (x): 

Y^x) := c )1 (x)n(x) + (?~ 1 9 At n(x) x n(x). (3.9) 

As a way of specifying the separation of variables: 

^(x)=V M (x)+X M (x), (3.10) 

the color field n(x) is required to be a covariant constant in the background field 
V(x): 

D$f\n{x) := d^n(x) + gY^x) x n(x) = 0. (3.11) 
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In fact, solving this equation for V M (x) leads to (3.9). On the other hand, the re- 
maining variable, i.e., the covariant potential X M (x) is required to be perpendicular 
to n(x): 

n(x)-X li (x)=0. (3.12) 

This decomposition (3.7) was once called the Cho-Faddeev-Niemi-Shabanov decom- 
positions [20-23] in the literatures. It is regarded as a non-linear change of variables 
(NLCV) for the original Yang-Mills field variables. 

For our purposes, it is a remarkable fact that the curvature tensor [V] obtained 
from the connection V M is parallel to n and its magnitude coincides exactly with 
the curvature tensor G^ u appearing in the Wilson loop operator (2.11) by way of the 
NAST: 

J^[V](2) :=d,Y u (x) - d u Y,(x) + gY^x) x Y u (x) 
:=n{x)G ^ u (x) = G^(x), 
G^ u (x) =d IM c v (x) - d u c^(x) - g~ l n{x) ■ (d^n(x) x d v n(x)). (3.13) 

Therefore, we have succeeded to separate the original variables 

K^x), n(x)) -> (c^x), X„(z), n(x)), (3.14) 

with the identification (3.8) such that only n(x) and c M (x) in the combined form V M 
are responsible for the Wilson loop average and that the remaining variable X M (x) is 
redundant for calculating the Wilson loop average. In other words, V M can be identi- 
fied with the "Abelian" part of =2^, suggesting the "Abelian" dominance in the Wilson 
loop average. This fact has been already pointed out in the paper [24]. However, this 
fact alone is not sufficient to guarantee the infrared Abelian dominance, since the the- 
ory has interactions between n, c M and X^. In order to confirm the infrared Abelian 
dominance, we must show that the variable X M (x) is actually irrelevant for calculating 
the Wilson loop average. This follows if these degrees of freedom decouple at least 
in the low-energy or long-distance region corresponding to the large Wilson loop. A 
possible mechanism is discussed in what follows. 

The SU(2) Yang-Mills Lagrangian density for the gluon field g/^, 

2ym\8*\ ■= -\{^A^\f = - \{d^u - + x sf v f, (3.15) 
is rewritten in terms of the new variables (n, c M , X M ) into 

& YM [n, c, X] = -i(G^ + gX^ x X,) 2 - i(I> M [V]X„ - A,[V]X, t ) 2 , (3.16) 

where we have used a fact that G^ and gX M x X„ are parallel to n, and this is also 
the case for the sum G^ + gX M x X y , while -D^fVjX^ — £)^[V]X M is orthogonal to n 
(which follows from the fact n ■ X M = 0). 

By collecting the terms in X M , the Lagrangian reads 

^ YM [n,c,X] = - X -Gl - \x$W$?X* - ^(gX, x X,) 2 , (3.17) 
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where 

■= ~ ^(D p [W]D p [W]) AB + 2ge ABC n c G, u + (D,,[V] A,[V])^. (3.18) 

In this derivation, we have used the relation: [.D^V], -D !/ [V]] AB = —ge ABC <S£ u . 

We return to the issue raised above: how to reduce the master Yang-Mills theory 
to the original Yang-Mills theory. If we treated the color field n(x) as the fundamental 
field in addition to the original gauge field £^(x), the resulting theory, say, the master 
Yang-Mills theory, had more independent degrees of freedom than those in the original 
Yang-Mills theory. In other words, the master Yang-Mills theory had the gauge 
symmetry G := SU(2) UJ x [SU (2) / U {1)]q larger than the original gauge symmetry G = 
SUityu- Here the latter symmetry [SU(2)/U(l)]g is carried by the color field n(x) 
while the former SU(2) W by the gauge field s^ p {x). In order to obtain the equivalent 
theory to the original Yang-Mills theory, we must impose necessary and sufficient 
numbers of constraints which eliminate the extra degrees of freedom and restrict the 
larger gauge symmetry to the SU{2) gauge symmetry. A suitable procedure was given 
in [25] by giving an explicit form of the constraint, which we called the new Maximal 
Abelian gauge (nMAG), although this naming is somewhat misleading. The nMAG 
is performed by minimizing the functional / d D x^g 2 lL 2 l with respect to the enlarged 
gauge transformations: 

= 6 U>9 J d D x^g 2 Xl = 8 ufi J d D x(D,[^]n) 2 . (3.19) 

This determines the color field n(x) as a functional of a given configuration of ^(x). 
The local gauge transformation of X 2 is calculated as [25] 

<J W ,^XJ = g-^D^n) ■ {5,K](wi - 0±) x n}, (3.20) 

and the average over the spacetime of (3.20) reads 

S u ,e J d D x^Xl = - J d D x(uj ± - ± ) ■ £>„[V]X„. (3.21) 

Hence, imposing (3.19) for arbitrary uj± ^ 6± yields a constraint in the differential 
form: 

X ■= D^[Y]X^ = 0. (3.22) 

This constraint yields two conditions, since n ■ £) M [V]X M = 0. Imposing nMAG to the 
master Yang-Mills theory breaks the enlarged G := S77(2) w x [SU(2)/U(l)]g gauge 
symmetry down to the diagonal SU(2) gauge symmetry: G' = SU{2)^i = SU(2)n, 
a subgroup of G (u> = := uj'). The respective new variable transforms under this 
gauge transformation, say the local gauge transformation II as follows [25]. 

Local gauge transformation II (the active or background gauge transformation): 

5' w n(x) =gn(x) x a/(x), (3.23a) 

S' ul c IM (x)=n(x)-d^'(x), (3.23b) 

c^X M (x) =gX p (x) x u'(x), (3.23c) 

^W(x) = £> M [V]u/(x), (3.23d) 

3X0*0 = ( 3 - 23e ) 



7 



Therefore, is invariant under the local gauge transformation II: 

O5(a0 = 0. (3.24) 

The gauge transformation for the field strength is calculated using this result. The 
curvature G^ := J^t/[V] is subject to the adjoint rotation 

S'J&^x) =gG^{x) x (3.25) 

Hence, the squared field strength has the SU(2) n invariance 

S'JZ^xf =0. (3.26) 

The inner product of G^ with n, i.e., the magnitude of G^ u , is also SU(2)n invariant: 

6' u (n(x) ■ G, u {x)) = SiG^x) = 0. (3.27) 

Thus the Abelian field strength G^ u is invariant under the SU(2) gauge transformation 
II, in sharp contrast to the original field strength which transforms in the adjoint 
representation. 

In the functional integral formulation, we must specify the integration measure 
for the new variables. We must take into account the constraints 5(n(x) ■ n(x) — 1) 
and S(n(x) ■ X M (x)) in the integration measure Pn(i)Dc (i (x)I'X (1 (x), 

dfi[n}V^(x) = Vn(x)5{n(x) ■ n(x) - l)Vc^x)VX^{x)5{n(x) • X^x)). (3.28) 

To avoid complications coming from constraints in performing the integration, we 
rewrite the integration measure in terms of the independent variables and calculate 
the Jacobian associated to this change of variables. For this purpose, we introduce the 
ortho-normal basis (ni(x), n 2 (x), n 3 (x)) = (ex(x), e 2 (x), n(x)), i.e., rij(x) ■ n k {x) = 
S jk , rijix) x n k (x) = e jki n £ (x), (j,k = 1,2,3), or equivalently 

e a {x) ■ e b {x) =8 ab , n(x) ■ e a {x) = 0, n(x) ■ n{x) = 1 
e a (x) x e b (x) =e ab n(x), n(x) x e a (x) = e ab e b (x), (a, 6 = 1,2). (3.29) 

The gauge transformation II for the basis vector rtj{x) is given by 

s L n j( x ) = 9nj(x) x u'(x) <^> ^e a (x) = ge a (x) x u>'(x), 5' u n{x) = gn(x) x u>'(x). 

(3.30) 

It is easy to show that the Jacobian is equal to one for the transformation 
from the original variables n A ,£/^ to the new variables n , c^X^ in this basis 
(ei(x), B2(x), n(x)) where (see Appendix A. 3 for the derivation): 

X M (x) = X;(x)e a (x) <-> Xfr) = X;(x)e^(x) 

or X«(x)=X ll (x)-e a (x)=X*(x)ei(x) (A = 1, 2, 3 : a = 1, 2), (3.31) 
so that the integration measure is written as 

dfi[n]V^{x) = Vn a (x)Vc^x)VX*(x). (3.32) 
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Thus, we have given a reformulation of Yang-Mills theory in term of new variables 
obtained by using the non-linear change of variables. In this reformulated Yang-Mills 
theory, the Wilson loop average is given by 



W(C) =Z^ M j Vn a (x)V Cfl (x)VX;(x)exp (^ YM [n, c, X]) W*{C) 
jVn a (x)Vc^{x)VX*(x) exp Us YM [n, c, X]) W^(C) 



(3 33) 

/ Vn a (x)Vc^x)VX^x) exp (iS yM [n, c, X]) 

where we have omitted the gauge-fixing term corresponding to nMAG (regarded as the 
condition for the partial gauge fixing from G to G 1 ) and the associated the Faddeev- 
Popov ghost and antighost term, see [26] for details. 

By making use of the gauge- invariant Abelian field strength G^ u , we can define 
the SU(2) gauge-invariant monopole current by 

k»(x) :=d v *GT(x) = (l/2)e^ pa d u G pa (x). (3.34) 

Then the magnetic charge is defined by the volume integral of k°: 

q m := / d 3 xk°. (3.35) 
Jv 

This is cast into the surface integral over the closed surface S as the boundary of the 
volume V, S — dV: 

q m = j> d 2 S/G oe = j> d 2 S e ^ k G jk = j da 3k G 3k . (3.36) 

It is easy to show that the gauge-invariant magnetic charge q m defined in this manner 
satisfies a charge quantization condition of the Dirac type: 

47T 

q m = — n (n = 0,±l,±2,.-.). (3.37) 
9 

Lattice formulations based on new variables of SU(2) Yang-Mills theory were for 
the first time constructed in [31,32] so that they reduce in the (naive) continuum 
limit to the reformulated SU(2) Yang- Mills theory written in terms of new variables 
obtained through NLCV from the original gauge field [26]. Two lattice formulations, 
non-compact [31] and compact [32], enable one to define the magnetic monopole in 
the gauge invariant way keeping color symmetry in Yang-Mills theory on a lattice 
without introducing fundamental scalar fields. This is a remarkable result, since 
the conventional approach of defining the magnetic monopole in Yang-Mills theory 
without fundamental scalar fields heavily relies on a specific choice of gauge fixing, the 
so-called the maximal Abelian gauge (MAG) [8]. The MAG breaks the local 577(2) 
gauge symmetry into the maximal torus group U(l) and simultaneously SU(2) color 
symmetry into U(l). Therefore, only U(l) local and global gauge symmetry remain 
if MAG is imposed. In the non-compact formulation [31], however, the magnetic 
charge resulting from the magnetic monopole defined in this way is not guaranteed 
to be integer-valued. The compact formulation [32] was constructed so that the 
magnetic monopole defined on a lattice by the similar way to (3.36) is integer-valued 
and satisfies the quantization condition (3.37). In fact, these features were directly 
confirmed by numerical simulations [32]. 
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4 Dynamical mass generation for the gluon fields 
through the vacuum condensation of mass 
dimension— two \X^\ 

We can introduce the gauge-invariant mass term which is invariant under the local 
SU(2) gauge transformation II as pointed out in [25]: 

^ m = ^M 2 X 2 . (4.1) 
This gauge-invariant mass term is rewritten in terms of the original variables stf^. 

^ m = l -M 2 x {^ - V,,) 2 = ^M 2 - Cp n + g-%n x nf 

' M 2 x {DJ^]n) 2 , (4.2) 



2<? 

under the understanding that the color field n is expressed in terms of the original 
gauge field ^ by solving the nMAG constraint. Therefore, (or c M and n) plays 
the similar role to the Stiickelberg field to recover the local gauge symmetry. Note 
that c M , n and X M are treated as independent variables after the non-linear change of 
variables and the mass term is a polynomial in the new variable X M , although they 
might be non-local and non-linear composite operators of the original variables stf^. 

The proposed mass term (4.1) or (4.2) for the gluon should be compared with the 
conventional gauge-invariant mass term of Kunimasa-Goto type [27]: 

^kg = Mhi{{^ - ig- l Ud^) 2 } = M 2 ti{{UD^]U ] ) 2 }, U(x) = e ~ ix(x)/v . 

(4.3) 

This mass term is non-polynomial in the Stiickelberg field x( x )- This fact makes the 
field theoretical treatment very difficult. 

We proceed to argue that there occurs a novel vacuum condensation of mass 
dimension-two for the field X M , i.e., (— X 2 ^ ^ 0, 2 and that the field X M acquires 
the mass dynamically through this condensation. A naive way to see this is to use 
the mean-field like argument or the Hartree-Fock approximation which leads to the 
gauge-invariant mass term for X M gluons: 



1 



2 



-X 2 \8 AB -(-XfXf 



X^ B = -MlX„ ■ X M , Ml = -a 2 /-X 2 



-X*H ■ ' M X = 



P 



(4.4) 



Now we make this idea more precise. First, we decompose the field X M into 
the background field X M and the quantum fluctuation field X M around it, i.e., X M = 
X M + X M . Then, we expand the Yang-Mills action for the Lagrangian density (3.17) in 



2 We adopt the Minkowski metric g^ v = diag(l, — 1, — 1, — 1). After the Wick rotation to the 
Euclidean region, the Minkowski metric tensor is replaced by —8^ = diag(— 1, — 1, — 1, — 1). 
Therefore, we have -X 2 -> (X^) 2 > 0. 
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X M around X M . The fact that the classical background field X M satisfies the equation 



of motion, ^xm^ c - : 



0, yields up to quadratic in X A 
1 



S YM [n, c, X] =5 YM [n, c, X] - ~X A X B K AB + (9(X 3 ), (4.5) 
where the quartic self- interaction term generates additional quadratic terms in X^: 
(^SymK c, X 



K AB 



1 

+ g 2 [g,J AB W - g„x A x B - 8 AB x c x° + 2X A x B - x B x A ] 



(4.6) 



with 3 W AB defined by (3.18). Hereafter, we neglect the quantum fluctuation parts of 
the other fields n(x) and c p (x) by identifying them entirely with low-energy slowly- 
varying modes: n(x) = n(x) and c^ix) = cf(x). This approximation will be improved 
by including also the high-energy modes for n(x) and c p (x) in a subsequent paper [28], 
although such an analysis has been tried in somewhat different context in [29]. 

By expanding the field X M (x) in the basis (ei(x), e2(x)) which is perpendicular to 
n(x), 

X M (x) = X;(x)e a (x), or X A (x) = X;(x)e A (x) (A = 1, 2, 3 : a = 1, 2), (4.7) 
the quadratic form is rewritten in terms of the independent fields X^(a = 1,2): 

\X A K AB X B = \x;Kf u X b v , Kf v := e A e B K AB . (4.8) 

Including the additional parts coming from quartic self-interactions among X M 
gluons is equivalent to modify the two-point gluon Green function (full gluon prop- 
agator) St^l by taking into account the tadpole contribution IT^ of the background 
field: 

K = HKr 1 + n ?- ( 4 - 9 ) 

{DfX 1 =Q%, K = V^K c p Xi, (4.10) 
where the X^-independent part Q^ u (x) is defined by 4 

-=\(W^(x) + W BA (x))e A (x)e B (x) = g, v R a \x) + 2gG fll/ (x)e ab , 
R a \ x ) :=5 ab [-d 2 + g 2 G p {x) 2 } + ge ab [ d p G p {x) + 2G p (x)d p ], (4.11) 



3 Kfif should be always understood in the quadratic form, \X^X® Kff . Therefore, \(W A „ + 
Wy A ) is equal to W^f under this understanding. 

4 In defining Q°fl, we have dropped the last term (^[Vjflj/JV]) coming from W^ B . If 
we introduce a gauge- fixing parameter a for nMAG, we have the gauge-fixing term of the 
form: / d D x±(D v \W\X v ) 2 = -/ d D x^X A {D^[¥]D"[¥}) AB X B which cancels the contribution 
^X A (D IJ '\V]D l '\y]) AB 'K^ for a = 1. Or. such a contribution vanishes by taking into account the 
nMAG condition _D"[V]X„ = 0, which corresponds to a = 0. Therefore, the following calculations 
should be understood to be performed exclusively in these cases. Moreover, it is worth remarking 
that exact satisfaction of nMAG condition (3.19) is realized only in the case a = 0. 
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and Vff%p is the four-point vertex for off-diagonal gluons [30] 

y abed 2r ab ed t , ac bdj , ad be j i j ._ I _ Wo 

(4.12) 

In deriving (4.11), we have used (3.29) to obtain 

D„[V]e a (aO = ^G At (x)e afe e b (x), (4.13) 

where G M is the SU(2) gauge-invariant Abelian gauge field defined by 

G>(x) = c M (s) + h^x), hpix) = g~ 1 d IJi ex{x) ■ e 2 (x) = -g~ l d^e 2 (x) ■ e^x). (4.14) 

These relations are shown to hold in Appendix A.l. Eq.(4.10) is nothing but the 
Schwinger-Dyson equation for the full gluon propagator in this approximation. This 
is a gap equation similar to that of four-fermion model of the Nambu-Jona-Lasinio 
type where the fermion mass is dynamically generated in a self-consistent way. 

Suppose that neither color symmetry nor Lorentz symmetry are spontaneously 
broken. Then we can not orient the component X® to a specific direction and the 
vacuum condensate is to be caused isotopically in color and Lorentz indices: 

K = V X (x^) = C 6 ^ C V <202S> = \a 2 9pJ ab (X?) • (4.15) 

Thus, we conclude that the existence of the vacuum condensation generates 
the mass term for the gluon field X M : 

~x;k;IxI = ~x;Qf u x b v + \m 2 x % • x", m\ : = ^ 2 (-x;) , (4.ie) 

where X 2 = X^X^ = X®X®. Then the X M gluon modes decouple in the low- 
energy (or long-distance) region below the mass scale Mx- Consequently, the infrared 
"Abelian" dominance for the large Wilson loop average follows immediately from the 
fact that the Wilson loop operator is written in terms of V M alone. 

The numerical simulations on a lattice [32] have demonstrated the infrared 'Abelian" 
dominance and magnetic monopole dominance in the string tension within our com- 
pact lattice formulation, although such phenomena were found for the first time in 
the MAG [5,6]. In fact, the string tension calculated from the magnetic part of 
the Wilson loop average according to (3.33), (3.3) and (3.4) in this formulation 
reproduces 90 ~ 95 % of the full string tension calculated from the original Wil- 
son loop average in the conventional lattice formulation. More numerical simula- 
tions [33] have shown that the remaining field X M defined on a lattice acquires the 
mass Mx — 1.2GeV which is obtained as the exponential decay rate of the two-point 
correlation function (X^(x) • X At (y)) YM measured on a lattice: (X At (x) • X At (y)) YM ~ 
\x — y\~ a exp(— Mx\x — y\). This value agrees with that of the off-diagonal gluon 
mass in the MAG [7]. On the other hand, the same analysis applied to the "Abelian" 
gluon leads to the result My = 0.6GeV. This is consistent with the "Abelian" 
dominance. 
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5 Effective potential and the vacuum condensation 



In order to study which type of the vacuum is realized in Yang-Mills theory, we need 
to calculate the effective potential. In particular, we pay attention to see whether 
such a dimension-two vacuum condensation occurs or not. First, we integrate out 
the fluctuation field X M in the functional integration. This is easily done for field 
with a trivial Jacobian for the change of variables: 



/ 2>X-<*p l-dx;K^X b v ] = (det^t)" V2 = exp 



-i— lndet^ 



Therefore, we are to calculate the effective potential 

V(X 2 ,G) = -G 2 + ^(X M x X,) 2 + ItrlnK 



ab 



-trln 

i 



ab 



4 -/«/ ■ 4 v-p ■ 2i 

where the quartic term in X M is decomposed into two gauge-invariant pieces: 

,2 



^(X M x X, 



fX M x X fc 



=^(VX") a -^(VX,) 2 . 



(5.1) 



(5.2) 



(5.3) 



Hereafter the underline for X^ will be omitted for simplifying the notation. The final 
term in (5.2) comes from the integration over the ghost and antighost fields which are 
necessary to implement the nMAG correctly according to the BRST method, see [26] 
for details. It should be remarked that this effective potential is SU(2) gauge invariant, 
since the Abelian gauge field G M is SU(2) gauge-invariant as well as the Abelian field 
strength G^, see Appendix A. 2. This is an advantage of our formulation, contrary 
to the conventional approach. 

For our purposes, we examine the effective potential as a function of two vacuum 
condensates X 2 , and G. 



1 



V(X\G) = l G% 



9 



2\2 



4 v ' 2i 

where K°^ v is shifted from Q ab u defined in (4.11) as 

3 



1 trln if J - -trlniT 6 , 



K 



ab 
fiu 



2^2 . 



;ab 



(5.4) 



(5.5) 



In calculating the effective potential, G^ v and X 2 are assumed to be constants uniform 
in space and time. 

In the limit of vanishing Abelian condensation G^ 
V(X 2 , 0) is written in the closed form: 



V(<f>,0) 



16 1 



In 



3 



0, the effective potential 
3 



2^2 



9 g 2 ^ (4vr) 2 \ ^ 
This is calculated as follows. The dimensional regularization yields 
lndet{[-<9 2 + 0]} = trln{[-<9 2 + 0]} 



(5.6) 



( 47r )D/2 

1 6 2 



r- 



r(-2 + e) 



(4vr; 



d.2-6 



2-e 



2 (4vr)' 



e" 1 + In 4tt - lE - In + - + 0(e) 



(5.7) 
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where e := (4 — D)/2 and 7s is the Euler constant. Then we obtain 



16 1 



2<f) 2 



9 g z (47T 



e" 1 + ln47r - lE - In + - + (9(e) 



We introduce renormalization constants for <fi and g as 



Then the effective potential reads 

16 1 n 2fh 

v(<P,o) = ^-<p 2 R z; 2 z, 

J 9r 



a/2 



2 7 



»fl> = Z g g R . 



(4tt) 



ln47r — 7s — In 



1/2 



From the observation, 



y r 



-2<f> R - -2 4> R (i - z g 2 z^), 

9r 9r 



and 



(4vr) 2 



(4vr) : 



ln47r — 7s — In 



W 2 



(4tt) 



/i 2 2 



(4tt) 



+ 



-In 



(4tt) 



M 2 2 



(4tt) 



e 1 + ln47r — 7s — In 



M 2 
7? 



e 1 + ln47r — 7s — In 



. 3 



1/2 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



these renormalization constants are chosen so that the counter term -^-(j) R (l — Z g 2 Z^) 

cancels the divergence in the bare effective potential, namely, to the lowest order of 
the coupling g R , 



16 1 M 7^7 



(47T) 



ln47r — 7s — In ■ 



M 2 
~7 



0. 



(5.13) 



where we have used Z^, Z g = 1 + 0(g 2 ). Thus we obtain the renormalized effective 
potential: 



hi 



9r_ 

M 2 



(5.14) 



which satisfies the renormalization condition: V"(cj)R = M 2 ,0) = (32/9)/g 2 . Indeed, 
the potential V{(t>R, 0) has a minimum at (f) R = <p R ^ away from the origin: 



<j>° R = M 2 e 1//2 



exp 



W)9 2 r(M) 



(5.15) 
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Figure 2: The effective potential V(4>, 0) where M = 1 and g = 10. 



See Fig. 2. This results shows that condensates indeed. In order to establish this 
condensation, this one-loop result should be improved. 

Operating the differential operator jjt,-^- to (5.13) and defining the /3-function P(gR) 



and the anomalous dimension 7^((7_r) as functions of g R by 



P(9r) = V 



dg R 



-gup- 



d\nZ n 



l<t>{9i 



1 d\nZ d 



<91n< 



dp dp 2 dfx 

we obtain a relationship between /3(gn) and 70(5^) up to C(g 3 ): 



/'- 



(5.16) 



R 



(47T) 



r<4 



(5.17) 



This is consistent with the asymptotic freedom [35], P(gn) = —j^^g%, with b = y 
provided that 



= 0h 



(4tt) 



;9r ~ P(9r) 



§ + b 2 
'9r- 



(47T) 



(5.18) 



In other words, the renormalized effective potential satisfies the renormlization group 
equation: 



M 



g 



P(9 



R 



d 

dg R 



10(9 r)4> 







R- 



>R 



RWR 



0. 



(5.19) 



Calculating the anomalous dimension 7^ in consistent with this relation will be given 
elsewhere [28]. 

This result should be compared with the effective potential V(a) calculated in 
[36] for the gluon-ghost mixed composite operator of mass dimension-two, ^A^A® + 
iaC a C a , which is shown [10] to be on-shell BRST and anti-BRST invariant in the 
modified MA gauge [37]. However, the vacuum energy E = V(ao) reached at the 
minimum depended strongly on the gauge-fixing parameter a of MA gauge. This 
indicates that the vacuum realized at cr does not corresponds to the true vacuum. 
Our effective potential V(<fr, G) is guaranteed to be gauge invariant by construction. 
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Therefore, we can choose any value of gauge fixing parameter for nMAG. In particular, 
the choice a = 1 simplifies the calculations so that the final term in cancels with 
a gauge-fixing term of nMAG. This choice is also preferred from a fact that the 
mixed composite operator for a = 1 reduces to the gauge-invariant part of the gluon 
composite operator j^A^A®, as demonstrated in [11]. The choice a = is the best for 
realizing the nMAG condition exactly. The true vacuum should be obtained at which 
the potential V(4>, G) takes the minimum. Thus the total effective potential could be 
gauge parameter independent, if the contribution from the Abelian part is included, 
even in the MA gauge. 

In the lattice formulation, the effective potentials for various composite operators 
of mass dimension-two [9-11] were calculated numerically on a lattice [34]. The re- 
sult shows that the numerically obtained effective potential for X 2 has a minimum 
away from the origin, suggesting the existence of non-vanishing vacuum condensates 
XjO 7^ 0- However, it should be remarked that they are unrenormalized quantities. 
We need to perform the non-perturbative renormalization to obtain the definite re- 
sult for the existence of such a dimension-two condensate. Therefore, it is not yet 
confirmed whether they survive in the continuum limit. 

6 Stability of the magnetic condensation 

Finally, we examine the contribution from the Abelian part G^ u = (E, H) . We 
distinguish two cases: (I) E • H ^ 0, (II) E ■ H = 0, characterized by two Lorentz 
invariants defined by 

T : = (E 2 - H 2 )/2 = —G^CT = (a 2 - b 2 )/2, (6.1) 
Q :=E H = —G^GT = ab, (6.2) 

where a = ^ T 2 + Q 2 +F, and b = ^V^ 2 + G 2 ■ 

(I) If E ■ H ^ 0, i.e., Q ^ 0, it is possible to transform to a Lorentz frame in 
which E and H are parallel or anti-parallel depending on the signature of E ■ H. 
We can choose the z axis as the direction of the vector without loss of generality: 
E = (0,0, E) and H = (0,0, H). The self-dual (or anti self-dual) case is a special 
case of (I): = *G^ (or G^ = -*G^), i.e., E = (0, 0, E) = H = (0,0, H) (or 
E = (0, 0, E) = -H= (0, 0, -H)). 

(II) If E ■ H = 0, i.e., Q = in a Lorentz frame, a = or b = 0; E and H are 
also perpendicular in any other Lorentz frame, since E ■ H is a Lorentz invariant. 
If E 2 > H 2 , i.e., T > 0, then the situation is as for the purely electric field with a 
vacuum instability, as is well known in QED. On the other hand, if H 2 > E 2 , i.e., 
T < 0, then the system behaves like the case of a purely magnetic field. If E 2 = H 2 , 
then we have a trivial case T = 0. 

The trace of the logarithm of a matrix is calculated, once all the eigenvalues of the 
matrix are known. In the pure magnetic case of (I), we can obtain the closed form of 
the effective potential, since the eigenvalues for R ab and Q ab u are exactly obtained to 
be capable of summing up all the contributions coming from all the eigenvalues. By 
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using the same method as in [38,39] and [40], we obtain 



V(<f>,H) 
1 , TT , 2 16 1 j2 



1 1 



(An) 2 4 



C(-i, 



— 1 4- r 1 4- r 

c(-i, — ^) + c(-i, ^— ) - c(-i, ^: 



+ 



c , (-i,^)+c(-i,4^)+c( 



ln^- + c 

1. 



gH' 

(6.3) 



where c := —1+7^ + In 2 with an Euler constant 7^ = 0.5772 ■ • ■ , and £(s, z) is called 
the generalized zeta function or Hurwitz (^-function defined by 



C(s, z) := £(n + > h ^0, -1, -2, 

n=0 

The generalized zeta function has an integral representation: 

1 



CM 



r z) J 



00 c~ z ^t s ~^~ 



dt- 



9£(s) > 1, 9fc(z) > 0, 



(6.4) 



(6.5) 



and its prime denotes the differentiation with respect to the first variable: £'(s, z) := 
j"C( s > 2 0- The zeta function is a special case of z = of the generalized zeta function: 
C(s) = C( s )0). The generalized zeta function can be analytically continued in the 
complex s plane to define an analytic function with a single simple pole at s — 1. 

In the limit r — > 0, i.e., neglecting the effect of vacuum condensation for off- 
diagonal gluons X, we recover the Nielsen-Olesen result [13]: 



v(o, h) = \h 2 



(4tt) 2 4 



where d :- 

c(-i,|) = 



c+ffC(-i,|) = 

-0.817409 • ■ ■ . Here we have used 



l+ 7B + ln2 + ffC / ( 



■1 -] 



-0.94556 



C(-i,|) = -i = C(-i,4) 



C(-i,|) + C(-i,^) 
C'(-i,^) + C'(-i,^) 



22 

=8[2C'(-l, |) - i|] = 16C(-1, |) - i47r. 



(6.6) 
with 

(6.7) 
(6.8) 
(6.9) 



This potential V(0, H) exhibits the instability due to the non-vanishing pure imag- 
inary part, j~g 2 H 2 . This is the so-called Nielsen-Olesen instability [13] to the 
magnetic vacuum of the Savvidy type [12]. 

The Nielsen-Olesen instability survives as long as r < 1. If r > 1, however, 
the pure imaginary part vanishes and the effective potential V((p, H) becomes real 
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number. This is because the four terms in front of In gH/fi 2 can be simplified by 
using 

(6.10) 



a-i,z) = - 1 -(z 2 -z+ 1 -), 



c(-i, m + c(-i, + c(-i, m - c(-i, |; 



22+9r 2 



and the primed 



2 / 1 'sV ^' 2 / 1 > 2 / ^ > 2' 

zeta function ('( — l,z) is real for z > 0. The last statement is checked by using the 
identity [41]: 

C'(-M) = C'(-l) " | ln(27r) - |(1 -z)+ £ dx\nT(x), (6.11) 

following from an integration of Binet's integral representation [42,43] of lnT(a;) and 
the Taylor expansion [42,44,45] of lnT(x) 



In T(x) = — In x — 7X + T"] 



■l) n 



n=2 



/i 



(6.12) 



where C'(-l) = C'(-l, °) = -0.1654 

The absolute minimum for V((j), H) exists in the region r > 1, and the value of 
r is determined as the point at which the potential takes the minimum. See the 
plot of V((j),H) in Fig. 3. Thus, the existence of dimension-two condensate (X 2 ) ^ 
guarantees as a by-product the stability of the magnetic condensation (H) ^ 
of the Savvidy type. In other words, the existence of dimension-two condensate 
(X 2 ) ^ shifts the gluon spectrum upward and eliminates the tachyonic mode causing 
the Nielsen-Olesen instability to recover the stability of the vacuum with magnetic 
condensation. 




Figure 3: The effective potential V{4>,H) vs. gH with a fixed value of r := 4>/{gH). (Left 
panel) r = 1.0, 1.5, 2.0, 3.0, 4.0, 5.0 from right to left. (Right panel) r = 10.0, 20.0. 



7 Conclusion and discussion 

In this paper, we have separated the original SU(2) gluon field variables ^ into 
and X M so that the variables V M are responsible for quark confinement and the 
remaining variables X M could decouple in the low-energy region. The former comes 
from a fact that a version of the non-Abelian Stokes theorem enables us to rewrite the 
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non-Abelian Wilson loop operator entirely in terms of the SU(2) invariant Abelian field 
strength defined from the variable V M . For the latter, we have argued that the 
gluon X M acquires the gauge-invariant mass dynamically through the non-vanishing 
vacuum condensation of mass dimension-two (X^ > 7^ 0. 



We have given a first analytical calculation of the effective potential of the com- 
posite gluon operator of mass dimension-two X^ demonstrating the occurrence of the 
condensation as realized at the minimum located away from the origin. See [34] for 
numerical simulations for the effective potential on a lattice. Note that the composite 
operator X^ is gauge-invariant and the resulting mass term for X M can be induced 
keeping the original SU(2) gauge invariance intact, contrary to the conventional wis- 
dom. Consequently, the decoupling of these degrees of freedom is characterized as a 
gauge-invariant low-energy phenomenon. In other words, this is a dynamical Abelian 
projection, suggesting the validity of the dual superconductor picture for quark con- 
finement. Thus the infrared "Abelian" dominance immediately follows in the gauge 
invariant manner 

A next issue to be investigated is to derive analytically the area law for the Wil- 
son loop average rewritten in term of new variables V M , since the area law with a 
string tension reproducing the full string tension was already shown numerically on a 
lattice [32], confirming the infrared "Abelian" dominance in the numerical way. The 
magnetic monopole is defined through V M (x) namely c M (x) and n(x) in the gauge- 
invariant way even in the Yang-Mills theory without introducing any scalar field as a 
fundamental field, just as the 'tHooft and Polyakov monopole in the Georgi-Glashow 
model. Therefore, the monopole dominance in the string tension can be in principle 
investigated, as confirmed by numerical simulations on a lattice [32]. 

Moreover, the existence of dimension-two condensate eliminates a tachyon mode 
causing the Nielsen-Olesen instability of the vacuum with magnetic condensation. 
Therefore, the restoration of the vacuum with magnetic condensation is obtained as a 
by-product of the above result. Furthermore, the existence of this condensate enables 
us to derive the Faddeev model describing glueballs as knot solitons as a low-energy 
effective theory of the Yang-Mills theory, as already pointed out in [38,46]. These 
are advantages of our reformulation of the Yang-Mills theory based on the non-linear 
change of variables. 
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A New variables in the ortho- normal basis (ei, e2, n) 

A.l Covariant derivative in the ortho-normal basis 

For Y p = c p n + g~ 1 d^n x n, we have 

Y p x ei = c M n x ei + g~ 1 (d fl n x n) x d 

= c M e 2 + g^fOVi ■ ei)n - d p n(n ■ d)] 

= c At e 2 + 5'~ 1 (9 M n-e 1 )n. (A.l) 

Then the covariant derivative of d in the background V M reads 

D fl [Y]e 1 := 9 M ei + x ei 

= c^d + gc p e 2 + (<9 M n ■ d)n 

= S'c^ea + <9 M d - (n ■ 9 M ei)n. (A.2) 

Now we define / M := d p e.\ — (^ - ^d) n - Then we can show easily that / M is orthogonal 
to d an d and hence / M is proportional to e 2 . Therefore / M is expressed as 

f t , = {f IJ ,-e 2 )e 2 = (e 2 -d p e l )e 2 =:gh p e 2 , = g~ 1 (e 2 ■ d^). (A.3) 

Thus we obtain 

D M [V]d = gc p e 2 + #/ve 2 = 5 , G M e 2 , G M = c M + V (A.4) 
In the similar way to the above, we can show that 

D p [Y]e 2 = -gG p e x . (A.5) 
If X M is written in terms of the orthonormal frame (d, d, 

X M = A;d+Ajd, (A.6) 

then we obtain 

D P [V]X, =D p [V](Xle 1 + Xf l e 2 ) 

=XlD p [Y] ei + X*D p [Y]e 2 + d p X\e x + «9 p A> 2 

=XlgG p e 2 - XlgG pei + c^d + <9 P ^e 2) (A.7) 

and 

- D p [Y]D p [Y]X p 

= - XlgG p D p [N\e 2 + X^G^Vjd 

- fl^^ej - D p [Y] [d p Xle 2 ] - d p [XlgG p ]e 2 + d^XfaGfa 
= - XlgG p D p [Y]e 2 + X*gG p D p [Y]ei 

- D p [Y][e x ]d p Xl - D p [Y][e 2 ]d p Xl - d p [d p X^e x - d p [d p X^e 2 - d p [XlgG p ]e 2 + d^XfaGfa 
=g 2 G p G p Xle l+ g 2 G p G p Xle 2 

- gG p n 2 d p X\ + gG p n x d p Xl - d p \d p X^e x - d p [d p X*]e 2 - d p [XlgG p ]e 2 + d p [XlgG p ]e x . 

(A.8) 
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Thus we obtain 

X p • [-D p \y\D p \V\]X„ 
=g 2 G p G p X p • X p + (-Xld 2 Xl - Xld 2 X 2 p ) 

- XlgG p d p Xl + XlgG p d p Xl - X 2 p d p (gG p XD + Xld p (gG p X 2 p ) 
=X«[-d 2 6 ab + g 2 G p G p 5 ab + 2ge ab G p d p + ge ab d p G p }X b 
=X+[-d 2 + g 2 G p G p + i(2gG p d p + gd p G p )]X; 

+ X;[-d 2 + £ 2 G P G P - i{2gG p d p + ^ P G P )]X M + 
=X+[-(5 p - igG p ) 2 ]X; + X;[-(«9 P + igG p ) 2 ]X+. (A.9) 



A. 2 New variables and conventional Abelian projection 

For the orthonormal basis (ei(x), e2(x), n(x)), we have 

X p =X 1 p e 1 + Xf l e 2 , X 1 p = e 1 -X p , X 2 = e 2 . X p , (A.10) 
B p =B\e x + 5 2 e 2 , B 1 ^ = g-%n ■ e 2 , B\ = -g~%n ■ e 1} (A.ll) 
/i p =g~ 1 d p e 1 ■ e 2 = -g~ l d p e 2 ■ e x . (A.12) 

For this basis, two gauge transformations are expressed as follows. 

Gauge transformation I: uj(x) = ui(x)ei(x) + u 2 (x)e 2 (x) + 9{x)n{x) 

<Lc„ =d p 6 + ig(X±cu* - Xju^), (A.13) 

S U X^ =[8 P + ig{c p + / ip )]^ ± =F z<?Xj0, (A.14) 

S U K =0, (A. 15) 
Gauge transformation II: uj'(x) = uj[(x)ei(x) + u' 2 (x)e 2 (x) + 9'(x)n(x) 

5^ =d p 6' + tg(B±Lu'T - Bju'^), (A.16) 
6„>X± =0, (A.17) 
Su'hp = - Su'ty, (A.18) 

where we have defined O := + i(9 2 ). 

In this basis, the nMAG reduces to the conventional MAG apparently: 

D P [Y]X^ = <^ d„X; - g(c„ + K)e ab X b p = 0, (A.19) 

provided that the Abelian part in the conventional Abelian projection is identified 

as 

a„ <- c M + h, : = G„, A a p <- X" (A.20) 

since D,jy\X v = e^X* - g[c p + /^)X 2 ] + e^X 2 + g{c^ + V)X*]. It should be 
remarked that a p is invariant under the SU(2) gauge transformation II: 

5^G,, = 0. (A.21) 

This is also the case for X?.: 

5^X; = 0. (A.22) 
This shows that this is different from the naive Abelian projection. 
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A. 3 Jacobian for the non-linear change of variables 

The Jacobian is calculated as follows. 5 By using the above bases, the field &/ A is 
decomposed as 



l„ + X, = c „n + (B; + X;)e a , < = c/ + (5; + lX. (A.23) 
We consider the change of 3D + 2 variables: 
«S n B ) -> (c, X, a , n c ) (a = 1, 2; A, S, C = 1, 2, 3; /x, !/ 



0. 



,£>-!). (A.24) 



Here n s and n c should be understood as denoting two independent degrees of free- 
dom obtained after solving the constraint n A n A = 1. However, if we choose specific 
components (directions), the color symmetry is apparently broken. Therefore, we 
keep this notation in the followings, keeping this convention in mind. 

The Jacobian is the determinant for the (3-D + 2) x (D + 2D + 2) matrix: 



dsrf., dn 

H 1 



B 



JdCydX^dn , 



J 



dc v dX£ !h£~ 
dn B dn B dn B 



(A.25) 



Since c u ,X^,n c are independent, we have 



dn 1 



0. 



dn 



B 



0, d ^ = 5 BC 



dc p dX* dn c 

Then the Jacobian reduces to the determinant for the 3D x (D + 2D) matrix: 



(A.26) 



J 



ds/A ds/f 

1 



Making use of (A.23), we have 



dc u 



and we conclude 
J 



n A e A 



dsrf A 
~dXj 



neie 2 



S e A 



\n ■ (e l x e 2y 



1. 



(A.27) 



(A.28) 



(A.29) 
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